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ABSTRACT 

In spite of a large number of global three-dimensional (3D) magneto-hydrodynamical (MHD) sim¬ 
ulations of accretion flows and jets being made recently, their astrophysical relevance for realistic 
situations is not well known. In order to examine to what extent the simulated MHD flows can ac¬ 
count for the observed spectral energy distribution (SED) of Sagittarius A* (Sgr A*), for the first time 
we calculate the emergent spectra from 3D MHD flows in a wide range of wavelengths (from radio to 
X-ray) by solving the 3D radiative transfer equations. We use the simulation data by Kato, Mineshige, 

& Shibata (2004) and perform Monte Carlo radiative transfer simulations, in which synchrotron emis¬ 
sion/absorption, free-free emission/absorption, and Compton/inverse Compton scattering are taken 
into account. We assume two temperature plasmas and calculate electron temperatures by solving 
the electron energy equation. Only thermal electrons are considered. 

It is found that the 3D MHD flow generally over-produces X-rays by means of bremsstrahlung 
radiation from the regions at large radii. A flatter density profile, p cx with a < 1, than that 
of the advection-dominated accretion flow (ADAF), p oc ^ is the main reason for this. If we 

restrict the size of the emission region to be as small as ^ 10 rg, where Ts is the Schwarzschild radius, 
the MHD model can reproduce the basic features of the observed SED of Sgr A* during its flaring 
state. Yet, the spectrum in the quiescent state remains to be understood. We discuss how to resolve 
this issue in the context of MHD flow models. Possibilities include modifications of the MHD flow 
structure either by the inclusion of radiative cooling and/or significant contributions by nonthermal 
electrons. It is also possible that the present spectral results may be influenced by particular initial 
conditions. We also calculate the time-dependent spectral changes, finding that the fluxes fluctuate 
in a wide range of the frequency and the flux at each wavelength does not always vary coherently. 

Subject headings: accretion: accretion disks — black hole physics — radiative transfer — Galaxy: 
center 


1. INTRODUCTION 

There is a long research history in the theoretical 
modeling of black-hole accretion flows. The standard- 
disk picture was first established by Shakura & Sunyaev 
(1973) after many attempts. In this model the gravita¬ 
tional energy is efficiently converted to radiation energy 
and is finally radiated away. Then, the disk was predicted 
to be luminous and relatively cold, exhibiting multi-color 
blackbody spectra (Mitsuda et al. 1984). The standard- 
disk model is widely accepted as a model for disks with 
moderately high accretion rates, M < Le/c^ with Le be¬ 
ing the Eddington luminosity and c being the light ve¬ 
locity (Esin, McClintock, & Narayan 1997). In fact, the 
calculated spectra based on the standard disk model fits 
well with the thermal component of the high-state of 
the black hole candidates (BHCs, Ebisawa 1999), includ¬ 
ing the specific temperature gradient (Mineshige et al. 
1994), and probably the big blue bump of the active 
galactic nuclei (AGNs, Shields 1978; Malkan 1983). 

On the other hand, when the mass-accretion rate is 
much less than the critical value of ^ Le/c^, the radia¬ 
tion loss is inefficient and thus the accretion flow becomes 
a radiatively inefficient accretion flow (RIAF). In such a 
situation, the thermal energy of the gas can be advected 
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inward towards the black hole without much being ra¬ 
diated away. This basic idea, now known as the notion 
of advection-dominated accretion flow (ADAF), was first 
presented by Ichimaru (1977) and has been investigated 
in quite a lot of detail since the 1990’s (Narayan & Yi 
1994, 1995a, b; Abramowicz et al. 1995; for reviews, see 
Narayan, Mahadevan, & Quataert 1998; Kato, Fukue, & 
Mineshige 1998). Since the ADAF can reproduce hard, 
power-law spectra in the X-ray range, as has been ob¬ 
served, this model was thought to be a good represen¬ 
tation of disks in low-luminosity AGNs (LLAGNs) and 
in the BHGs during their low-hard state. Remarkably, 
the ADAF model can nicely fit the emergent spectrum 
of Sagittarius A* (Sgr A*) and other black-hole objects 
(Narayan, Yi, & Mahadevan 1995; Manmoto, Mineshige, 
& Kusunose 1997; Narayan et al. 1998; Manmoto 2000; 
Ozel, Psaltis, & Narayan 2000; Oka & Manmoto 2003). 
(We should keep in mind, however, that most of obser¬ 
vational data of Sgr A* only gives upper limits except at 
radio and X-ray wavelengths.) 

We should be aware that there are a number of seri¬ 
ous problems inherent to the ADAF formulation. For 
example, the ADAF model cannot properly treat three- 
dimensional motion, because the ADAF model is a one¬ 
dimensional model (in the sense that only the radial 
structure is solved), although both of the simple argu¬ 
ments and hydrodynamical simulations have shown that 
the occurrence of convections and/or outflow seems to 
be a natural consequence in the RIAF (Narayan & Yi 
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1994; Igumenshchev & Abramowicz 2000). The magnetic 
fields are treated as a single parameter, the plasma-/?, to 
predict synchrotron emissivity, which is too simple and 
problematic to describe the dynamics of magnetic fields 
(see a comprehensive discussion in Narayan 2002). Yuan, 
Quataert, & Narayan (2003) recently succeeded in repro¬ 
ducing the spectrum of Sgr A* by introducing nonther- 
mal electron components, but they still treat the mag¬ 
netic fields in a simplihed way. Recent MHD simulations 
have commonly shown the important role of the magnetic 
fields in the flow dynamics. The MHD flow is intrinsi¬ 
cally time varying and exhibits fractal structure (e.g., 
Kawaguchi et al. 2000; Machida & Matsumoto 2003, 
hereafter refereed to as MM03). To summarize, the lack 
of three-dimensional (3D) motion and the over-simplified 
treatment of magnetic fields lead to self-inconsistencies 
in the ADAF solution. 

The global 3D MHD simulations of RIAFs were hrst 
made by Matsumoto (1999), and have been extensively 
performed recently by several groups (see Mineshige & 
Makishima 2004 for a compilation of recent works). It 
has been revealed by the 3D MHD simulations that the 
flow pattern is considerably complicated and differs sig¬ 
nificantly from that of the ADAF model (Hawley 2000; 
Machida, Hayashi, & Matsumoto 2000; Hawley & Krolik 
2001; Machida, Matsumoto & Mineshige 2001, hereafter 
refereed to as MMMOl; Hawley & Balbus 2002, hereafter 
HB02; MM03; Igumenshchev, Narayan, & Abramowicz 
2003, hereafter INA03). Stone & Pringle (2001, here¬ 
after SPOl) established that a complex flow pattern is 
produced due to a magneto-rotational instability (MRI, 
see also Hawley, Balbus, & Stone 2001; Hawley 2001; 
Balbus 2003 for a review). Since the dynamics of mag¬ 
netic fields is fully solved, the 3D MHD models seem 
to be more suitable for describing the RIAF than other 
models without appropriate treatment of magnetic helds. 
Thus, the MHD models are expected to fit the obser¬ 
vations. Surprisingly, however, such a critical test has 
not been well investigated. Some papers discussing the 
observational appearance of the simulated MHD flow, 
HB02 for example, have shown the distribution of the 
synchrotron emissivity based on their MHD model. Mi¬ 
neshige et al. (2002) reported a preliminary examination, 
claiming that the radial density profile is significantly 
flatter than that of ADAF, leading to over-production 
of X-rays, but they calculated the spectra in a simpli¬ 
fied way. Their results should be confirmed by full 3D 
radiative transfer calculations. More recently, Goldston, 
Quataert, and Igumenshchev (2005) calculated the spec¬ 
tra based on the simulated MHD flow by INA03, but 
only in the radio band. To our best knowledge nobody 
has yet calculated the emergent spectra in a wide range 
of wavelengths (from radio to X-ray). This prompted 
us to calculate, for the first time, the emergent spectra 
based on simulated 3D MHD flows by performing 3D 
Monte Carlo radiative transfer simulations and directly 
comparing them with the observed data of Sgr A*. 

The plan of this paper is as follows: We present our 
model and the method of Monte Carlo radiative transfer 
in §2. The results will be displayed with the observational 
data in §3. We will demonstrate a serious discrepancy 
between them. A discussion of the loopholes, which cause 
the discrepancies, and how to remove them will be given 
in §4. The final section is devoted to conclusions. 


2. MHD MODEL AND METHOD OE RADIATIVE 
TRANSFER CALCULATIONS 

2.1. Overview of Adopted MHD simulations 

Our spectral calculations are based on 3D MHD sim¬ 
ulations by Kato, Mineshige, & Shibata (2004: hereafter 
refereed to as KMS04). They investigated the evolu¬ 
tion of a torus threaded by weak localized poloidal mag¬ 
netic fields. The Overall evolution of 3D MHD accre¬ 
tion flows is divided into two distinct phases (see Fig¬ 
ure 1 in KMS04). In the first phase (t<1800), toroidal 
magnetic fields are generated by differential rotation 
and are accumulated in the central region, driving a 
magnetic tower jet, where t is the time normalized by 
Ts/c ~ {M/Mq)s with M being the black-hole mass 

and Ts being the Schwarzschild radius. The jet is, how¬ 
ever, a transient phenomenon. Eventually, {t > 1800), 
the jet ceases and a quasi-steady, geometrically thick 
density distribution with complex held conhguration is 
produced. 

In the quasi-steady state, the radial density prohle is 
p ex r in the inner part (r < 20rs), while p ex r~^ in the 
outer part (r > 20rs) (see Figure 4 in KMS04). That is, 
there is a broad peak at r ~ 20rs in the density prohle, 
which might be a remnant of the initial torus (discussed 
later). Note that this density prohle differs from the 
results of some other 3D MHD simulations (see §4 for 
discussion). Also note that the density distribution of 
the ADAF model is signihcantly steeper; p oc with 
a = I-3/2. 

Throughout the present study, we set the black-hole 
mass to be M = 2.6 x IO^Mq (Schodel et al. 2002; Chez 
et al. 2003). Thus, the normalized time corresponds to 
t/i =26 s. We basically use the data of quasi-steady state 
at t = 2210 {t = 5.7 x lO'^s), when no obvious outhow is 
observed (see Figure 1 in KMS04), except at §3.3, where 
we will examine spectral variations. 

2.2. Calculations of Physical Quantities 

In the present study, we use Cartesian coordinates, 
{x,y,z), where the black hole is located at the origin of 
the coordinate axes, the z-axis is set to be the rotation 
axis of the accretion flow, and the x-y plane corresponds 
to the equatorial plane. We employ Cartesian grids with 
numbers Ny, Afz)=(100,100,197) of cells. We as¬ 
sume that radiation is generated within the cylindrical 
region with radius, i?; (x^ +2/^)^^^ < R and z < 100 Cg. 
We consider the cases oi R = 10 Xg and 30 rg. The size 
of the calculating box is 2X x 2Y x 2Z, where we set 
{X,Y,Z) = (R,i?,100 rg). 

We take 3D data of density, magnetic Helds, and pro¬ 
ton temperature distributions in the accretion flows from 
the 3D MHD simulations (KMS04). Since the MHD sim¬ 
ulation only give the normalized density, p, with the nor¬ 
malization, pq, and the normalized field strength, B, we 
have one free parameter, po, to determine absolute val¬ 
ues of density and magnetic fields; that is, p = pop and 
B = {poc^)^^^B, where B is the field strength. The pro¬ 
ton temperature does not depend on the density parame¬ 
ter and is given by the MHD simulation as /k)cs‘^, 

where p is the mean molecular weight (=0.5), Wp is the 
proton mass, k is the Boltzmann constant, and cj, is the 
normalized sound velocity obtained by the simulation. 
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Although KMS04 assumed a one-temperature plasma, 
here we adopt two-temperature assumptions. Assuming 
that the electrons have a Maxwellian distribution, we 
evaluate the electron temperature, Te, through the en¬ 
ergy balance of the electrons between Coulomb collisions 
with ions and radiative cooling, 

rZ pR r 


/ / Xie^Trrdrdz 


Ludv. 


( 1 ) 


J-zJo j 

Here, Aje is the energy transfer rate from ions to elec¬ 
trons (Stepney & Guilbert 1983), is the luminosity 
at frequency, respectively. In the present study, we 
suppose, for simplicity, the electron temperature to be 
a function of only the radius, r, and to be independent 
of the altitude, z. This simplification does not affect 
the results too much, since the emission from regions 
around the equatorial plane, where the density is at a 
maximum, is dominant at each radius. In addition, we 
solve for the electron temperature by dividing the cal¬ 
culating box into two parts; the inner (r < 10|rs) and 
outer parts (r = 10 — 30 Tg)- The fractional luminosity, 
Ly, is obtained by Monte Carlo simulations (see next 
subsection) for a given Te, which does not always satisfy 
equation O- We thus iteratively calculate the electron 
temperature and the emergent spectrum so as to meet 
the condition O- 


2.3. Monte Carlo Radiative Transfer Calculations 


The method of the Monte Carlo simulation is based 
on Pozdnyakov, Sobol, & Sunyaev (1977). In this 
study, synchrotron emission/absorption, free-free emis¬ 
sion/absorption, and Compton/inverse Compton scatter¬ 
ing are taken into account. Some of the photons emit¬ 
ted by synchrotron and free-free emission pass through 
the calculating box without being scattered or absorbed, 
while some are scattered and/or absorbed in the accre¬ 
tion flow, depending on the mean free path at the pho¬ 
tons emitting position. In order to efficiently calculate 
the emergent spectra, we introduce a weight, w, as de¬ 
scribed by Pozdnyakov, Sobol, & Sunyaev (1977). At the 
beginning, the weight, wq, is set equal to unity for each 
emitted photon and then we calculate the escape proba¬ 
bility, Po- The escape probability of a photon after the 
i-th scattering (for i > 1), P^, is evaluated as 


Pi = exp 


/ p{xi,yi,Zi) \ 

\ TUp j 


crK^{xi, yi,Zi) 





( 2 ) 


where {xQ,yo,Zif) corresponds to the point where a pho¬ 
ton was generated, {xi, yi, Zi) is the point where a photon 
is subject to the Tth scattering, Wp is the proton mass, 
(Jkn is the Klein-Nishina cross section (Rybicki & Light- 
man 1979), Kabs is the absorption coefficient, and I is the 
distance from the point {xi, yi,Zi) to the boundary of the 
calculating box, respectively. Assuming local thermody¬ 
namic equilibrium (LTE), we can write the absorption 
coefficient as 


K 


abs 

v 


47rP;, 


( 3 ) 


where and e® are the synchrotron and free-free emis- 
sivity, respectively (Pacholczyk 1970; Stepney & Guil¬ 
bert 1983) and is the Planck function. The quantity 


of wqPq represents the transmitted portion of photons 
and is recorded to calculate the penetrated spectrum or 
reprocessed photons according to the escape direction of 
the photon. A fraction of Aq of the remaining portion, 
wi = wo(l — Po)Ao, undergoes at least one scattering, 
while u>o(l — Po)(l — Aq) describes the absorbed portion, 
where Ai {i > 0) is the scattering albedo. The transmit¬ 
ted portion of photons after the Tth scattering, WiPi, is 
recorded to calculate the transmitted spectrum, and the 
remaining portion, Wi{l — Pi)Ai, undergoes the (i-l-l)-th 
scattering. This calculation is continued until the weight 
Wi becomes sufficiently small {wi 1). The whole pro¬ 
cess is simulated by the Monte Carlo method. Finally, 
we suppose the region within the Schwarzschild radius, 
{x"^ Py"^ + z^Y!"^ < Ts, to be vacuum. Gravitational lens- 
ing as well as photon redshifts are not considered, since 
the black hole is much smaller than the calculating box 
in our study. 

3. RESULTS: SPECTRAL FEATURES OF MHD FLOW 
3.1. Basic features 

We first display the representative spectra of the MHD 
flow in Figure 1 (see thick solid curves) together with the 
observed data of Sgr A*. The adopted parameters are 
the radius of the emitting region, R — 30 rg, the elapsed 
time at t = 5.7 x lO^s, and the density normalization, 
Po = 1.6 X 10“^'^gcm“^ (upper) and 5.6 x 10“^®gcm“^ 
(lower). It might be noted that the emission is predomi¬ 
nantly from the equatorial plane and other parts at large 
altitudes (|z| 10 rg) do not contribute very much to 

the emergent spectra, although we sum up all the contri- 



Fig. 1.— The emergent spectra of the MHD accretion flow and 
the observed data of Sgr A*. The thick solid curves are the resul¬ 
tant spectra for R = 30 Vs. Here, the adopted density parameters 
are po = 1.6 x 10“^‘^g cm“^ (upper) and 5.6 x 10~^®g cm~^ (lower), 
respectively. The thin solid curve indicates the spectrum, where 
the magnetic fields are set to be ten times stronger than the spec¬ 
trum represented by the lower thick-solid curve. The spectra from 
the innermost part within 10 Vs are plotted by the dashed curves 
[Pq = 1.6x 10~^^gcm“^ (upper dashed), 2.0 x 10~^^gcm“^ (lower 
dashed)]. The filled circles and the lines with arrows indicate the 
data and upper limits by radio and IR observations. For more 
detailed information regarding the data, please refer to Narayan 
et al. (1998). X-ray observations of the fiaring and the quiescent 
state are shown by the two ‘bowties’ (Baganoff et al. 2001, 2003). 
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TABLE 1 

Calculated electron temperature. 


PO 

Te(r < 10 ra) Te(10 Ts < r < 30 rfl 

1.6 X 

5.56 X lO^'K 

3.45 X lO^K 

5.6 X 

4.93 X 10®K 

3.58 X 10®K 

2.0 X 

5.20 X lO'^K 

— 

3.0 X 10-l®gcm-3(Our 2D) 

4.21 X lO'^K 

3.65 X 10®K 


Our 2D: our 2D MHD simulations 

buttons from the vertically elongated cylinder between 
Z = ±100 Tg. This is because the density rapidly de¬ 
creases as \z\ increases. The calculated electron tempera¬ 
tures are listed in Tabled As expected, electron temper¬ 
atures are insensitive to radius, but there is a slight ten¬ 
dency that Te decreases with an increase of radius. Simi¬ 
larly, relatively flat electron temperature profiles were ob¬ 
tained in two-temperature, hot accretion flows (Narayan 
& Yi 1995b; Nakamura et al. 1996; Manmoto, Mineshige, 
& Kusunose 1997) and evaporated disk-corona flows (Liu 
et al. 2002). 

It is apparent that the resultant spectra look similar 
to those of the ADAF model. The lower-energy peak in 
the radio band is a synchrotron peak created because of 
significant self absorption of synchrotron emission. The 
IR emission at around log v = lA is, due to inverse Comp¬ 
ton scattering of synchrotron photons. Since we assume 
thermal electrons only, the spectral slope at logj^ < 11 is 
Lj/ oc corresponding to that of Rayleigh-Jeans emis¬ 
sion. The electron temperatures do not vary much, even 
if we change poi and, hence, the radio parts (which de¬ 
pends primarily on the electron temperature) are roughly 
identical among different models, as shown in Figure 1. 
In contrast, the flux at logz/ > 12 and the frequency 
of the lower-energy peak (in the radio band) are both 
sensitive to po- 

Here, we stress that the magnetic field strength affects 
the lower-energy peak in the radio band and the IR flux, 
whereas the X-ray flux and the slope at log v < 11 do not 
depend on the magnetic field strength itself. To check 
what alters if the magnetic field strengths were under or 
over-estimated in the simulation, we calculate the spec¬ 
tra with artificially strengthened magnetic fields, keeping 
the same density profile. The thin solid curve in Figure 1 
indicates the spectrum, here the magnetic field strength 
is set to be ten times stronger than the original values, 
so as to be compared with the original (lower thick-solid 
curve). It is found that the X-ray flux and the slope 
at logz^ < 11 do not change, but the radio peak lumi¬ 
nosity and IR flux both increases. This is because the 
synchrotron emission contributes to the radio flux and 
inverse Compton scattering of the synchrotron photons 
is dominant in the IR band. 

A big distinction between the spectra of ADAF and 
those of simulated MHD flows is found in the X-ray 
bands; namely, inverse Compton scattering is dominated 
in the former, while the thermal bremsstrahlung is the 
dominant mechanism in the latter. This reflects a flat¬ 
ter density profile in the MHD flow than for the ADAF, 
since a flatter density profile means that more material 
is present at large radii than the for case with a steep 
density profile, thereby producing more bremsstrahlung 
photons (discussed later). 


We confirmed that the resultant spectra do not signif¬ 
icantly change even if we employ the density and mag¬ 
netic fields averaged over the azimuth in the simulated 
data of KMS04. This implies that 3D effects are not es¬ 
sential for the study of the spectra. This is because the 
flow pattern is almost axisymmetric in the simulations of 
KMS04, who calculated the evolution of a torus threaded 
by weak poloidal magnetic fields. In the case that the 
perturbations in the azimuthal direction increase, and 
non-axisymmetric structures form, any 3D aspects would 
clearly appear in the spectra. 

3.2. Fitting to the flaring-state spectrum 

A striking fact is that we cannot fit both the radio and 
X-ray data simultaneously with the current MHD flow 
model. The MHD flow model can reproduce only a part 
of the observations. For example, it can fit the observed 
radio peak, if we assign po = 1-6 x 10“^'^gcm“^, but 
its flux largely exceeds the observed X-ray data and the 
upper limits in the IR band. This X-ray excess is caused 
by the strong free-free emission from the outer part of the 
flow (r > 10 Tg), since the emissivity of free-free emission, 

which is dominant in X-rays, is e® oc oc for a 

density profile of p oc r and the Tg profile of Tg oc r“^. 
The entire luminosity is s^d^r oc r^'^dr. 

For comparison, we plotted the spectra for the emis¬ 
sion from only the inner region within i? = 10 Tg in Fig¬ 
ure 1 (see the upper dashed curve) for the same density 
parameter, po = 1.6 x lO^^'^gcm”^. The electron tem¬ 
perature is 5.56 X 10®K (see Tabled). As shown in this 
figure, huge X-ray excesses disappear (cf. the case with 
i? = 30 Tg) and the model is successful in reproducing 
radio observations at around logzz 11 — 12 and is con¬ 
sistent with X-ray observations during the flaring state 
(see the upper dashed curve). 

The discrepancy between the upper thick-solid and up¬ 
per dashed curves represents the contribution from the 
outer parts (r > lOrg). This clearly demonstrates that a 
huge X-ray excess is generated at the outer part of the 
flow and we need to somehow remove these contributions 
to fit the date. 

Alternatively, we can reduce the X-ray flux so as to 
fit the data by employing a lower density parameter, 
Po = 5.6 X 10“^®gcm“^ (see the lower thick-solid curve in 
Figure 1), but for such a case, radio peak flux decreases in 
accordance with the reduced density and, hence, is short 
of the radio luminosity. We thus see that for this partic¬ 
ular density profile, it is totally impossible to reproduce 
the observations regardless of po (see also Mineshige et 
al. 2002). 

This situation resembles that of the CDAF 
(convection-dominated accretion flow; see Ball, Narayan, 
& QuataertBall 2001), since it has a flatter density 
profile like for the MHD flow (MMMOl). For p oc 
and T cx r~^, we have e^d^r oc r^'^dr, indicating 
significant free-free emission from the outer parts. 

Here, we stress that the emergent spectra starts to 
resemble to the observed data if the magnetic field 
strengths are systematically under-estimated in the 3D 
MHD simulation, that is, the actual magnetic fields are 
stronger than in the simulated data. As was already men¬ 
tioned in the previous subsection, the radio peak lumi¬ 
nosity and IR flux both increases with an increase of the 
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magnetic field strength. Therefore, the emergent spectra 
are consistent with the observations if the actual mag¬ 
netic fields are ten times stronger than the simulated 
data (see the thin solid curve). 

Strictly speaking, the calculated spectrum also devi¬ 
ates from the observed data at lower frequencies, log v < 
10. This inconsistency would be resolved if we consider 
synchrotron emission from nonthermal electrons as was 
first claimed by Mahadevan (1998). By a careful fit¬ 
ting to the up-to-dated spectrum. Yuan, Quataert, & 
Narayan (2003) concluded that the low-frequency radio 
flux can be explained via the emission of nonthermal elec¬ 
trons, which possess a small fraction of the electron ther¬ 
mal energy. 

3.3. Case of the quiescent spectrum 

What about the case of the quiescent spectrum? Be¬ 
fore attempting any fitting to the quiescent spectrum, 
we need to remark on the observational constraints of 
the size of the emission region. By the Chandra X-ray 
observations, Baganoff et al. (2001, 2003) have claimed 
that the X-ray emission of Sgr A* is extended on scales 
of about 0.04 pc in the quiescent state. However, Tan & 
Draine (2004) suggested that about half of the extended 
X-ray emission is due to dust scattering from an unre¬ 
solved source. Therefore, we regard the observed X-ray 
data in the quiescent state to give an upper limit for the 
X-ray luminosity of the accretion flow around the black 
hole. 

To fit the X-ray flux in the quiescent state, we need to 
adopt a small density parameter, po = 2.0 x 10“^®gcm“^ 
(see the lower dashed curves), but, then, the flow can not 
reproduce a radio peak. In addition to this, the X-ray 
spectral slope also differs. The X-ray excess decreases 
with a decrease in size of the emission region, since the 
X-ray emission is dominated by the outer part, as we 
already mentioned in the previous subsection. We can 
adjust the X-ray luminosity to the observed value in the 
quiescent state by setting a much smaller emitting region, 
R < 10 Tg. However, the X-ray spectral slope is still much 
steeper than the observed one. 

To summarize, we can reproduce the observed SED 
during the flaring state if the emitting region of the flow 
is restricted to be relatively small (r < 10 rg), but cannot 
account for the quiescent-state spectrum for any choice 
of parameters. Note that a compact emission region is 
needed to account for short-term variability, since oth¬ 
erwise the variability timescale will be much longer in 
accordance with the timescales in the outer zone (Mi- 
neshige et al. 2002). 

3.4. Time variation 

MHD simulations show that MHD flows are highly 
time-variable. It is thus interesting to see how the re¬ 
sultant spectra change with time and check if the calcu¬ 
lated trend fits the observational ones. We plot the time 
variations of the spectra in Figure 2. Here, we again 
set i? = 10 Tg and po = 1-6 x 10~^‘*gcm“^. The calcu¬ 
lated epochs are divided into two phases; the jet phase 
(t <5 X lO^s) and quasi-steady phase (t>5 x lO^s). The 
fluxes fluctuate in a wide range of frequencies as shown 
in the upper panel, although the electron temperature 
is kept nearly constant in the present study (see lower 
panel). 
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Fig. 2.— Spectral variations. Upper panel: the resultant spectra 
at t = 3.1 X lO^s (solid), t = 4.0 x lO^s (dotted), t = 4.7 x lO^s 
(dashed), t = 5.7 x lO^s (dot-dashed), and t = 10.1 x lO^s (long- 
dashed), respectively. Here, we set R = 10 rs and po = 1.6 x 
10“i^gcm~^. The lower panel: time variations of radio, IR, and X- 
ray luminosities as well as the electron temperature. The electron 
temperature is nearly constant. The radio and IR luminosities 
monotonously increase until the jet is produced (t < 5 x lO^s) while 
X-rays fluctuate. Radio luminosity is nearly constant at t > 5 x 
lO^s, but the flux slightly fluctuates at IR and X-ray wavelengths 
and their amplitudes vary a few hundred percent. 

As shown in this figure, the radio and IR luminosity 
increases until the jet is generated (t<5 x lO'^s), since 
magnetic fields are steadily being amplified in this phase. 
In the quasi-steady phase (t>5 x lO^s), in contrast, 
the radio luminosity is kept nearly constant, whereas IR 
and X-ray luminosities fluctuate and their variation am¬ 
plitudes amount to a few hundred percent. These re¬ 
sults are roughly consistent with the observed trends of 
multi-wavelength variability properties reported by Ul¬ 
rich, Maraschi, & Urry (1997), who show that the largest 
amplitude variations are in the soft X-rays and that the 
optical-UV variation amplitude decreases systematically 
with increasing wavelength in the LLAGNs. 

Interestingly, the calculated flux at each wavelength 
does not always vary coherently. We found that the X- 
ray luminosity decreases until t ~ 4 x lO'^s, while the 
radio and IR luminosities increase. Moreover, only the 
radio luminosity does not vary so much between t = 4.7 x 
lO'^s and t = 5.7 x lO'^s in spite of the significant decrease 
in the IR and X-ray luminosities (see lower panel). A 
comparison of such trends with the observations is left 
for future work. 

The X-ray time-variation predicted by our results is too 
small to account for the observed large-amplitude varia¬ 
tions (X-ray burst) reported by Baganoff et al. (2003), 
i.e., the X-ray luminosity in the flaring state is about 45 
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times larger than that in the quiescent state. Such an X- 
ray flaring event might be produced by the nonthermal 
electrons, whereas we take only thermal electrons into 
account in the present study. 

4. DISCUSSION: THE IMPLICATIONS OF OUR RESULTS 

We have calculated the emergent spectra based on the 
MHD simulation by KMS04 and demonstrated that the 
MHD flow model cannot account for both of the radio 
and X-ray observations in the flaring state simultane¬ 
ously, unless we restrict the emission region to be com¬ 
pact. We have also shown that it is problematic for re¬ 
producing the spectrum in the quiescent state. In this 
section, we discuss what these results imply. 

4.1. flaring state 

We first focus our discussion on the case of the flaring 
state. As was already mentioned in §3.2, the relatively 
compact emission region (r < 10 Cg) is required for the 
MHD flow model to reproduce the observed spectrum in 
the flaring state of Sgr A*. There are several possibilities 
which meet this requirement. 

4.1.1. Modifying density profiles 

A compact emission region would be realized, if the 
density profile at large radii would be much steeper than 
that obtained by KMS04 or if the broad density peak 
found in KMS04 is somehow removed. Since KMS04 pro¬ 
vides only one specific example, it is necessary to check 
the density profiles of other MHD simulations. In Table|21 
we summarize the density profiles of some representative 
simulations with brief descriptions of their calculations. 

SPOl, HB02, and KMS04 adopted similar initial con¬ 
ditions; namely, they started calculations with a torus 
threaded by weak localized poloidal magnetic fields and 
assumed no further mass input. They all investigated the 
evolution of the torus under the pseudo-Newtonian po¬ 
tential. INA03 and Proga & Begelman (2003: hereafter 
refereed to as PB03) set a different situation; INA03 con¬ 
tinuously injected magnetized matter into the computa¬ 
tional domain from near the outer boundary, and PB03 
started the simulations with Bondi accretion flow, whose 
angular momentum is zero except for the outer part of 

TABLE 2 

List of MHD models 


ref. 

'I' 

I.C./B.C. 

jet 

p-profile 

p-peak 

SPOl 

PN 

Sp-torus 

yes 

r^-^(< 10 rg), r“^(> 10 rg) 

10 rg 

MMMOl 

N 

B^-torus 

no 

r-l/2 

— 

HB02 

PN 

Sp-torus 

yes 

r®-)-torus(r < 10 rg) 

5 rg 

INA03 

PN 

injection 

yes 

10 rg), r“l(> 10 rg) 

3 rg 

MM03 

PN 

S^-torus 

no 


3 rg 

PB03 

PN 

Bondi flow 

yes 

r°-2(< 7 rg), r-l(> 7 rg) 

7rg 

KMS04 

PN 

Bp-torus 

yes 

rl(< 20 rg), r“l(> 20 rg) 

20 rg 

Our 2D 

PN 

Bp-torus 

yes 

rl(< 7 rg), r-l(> 7 rg) 

7rg 


SPOl: Stone & Pringle (2001); MMMOl: Machida, Matsumoto 
h Mineshige (2001); HB02: Hawley & Balbus (2002); INA03: Igu- 
menshchev, Narayan, & Abramowicz (2003); MM03: Machida & 
Matsumoto (2003); Our 2D: our 2D MHD simulations; 'L: gravita¬ 
tional potential; N: Newtonian; PN: pseudo-Newtonian; I.C./B.C.: 
initial or boundary condition; or Rp-torus: A torus threaded 
by weak localized toroidal or poloidal magnetic fields is initially as¬ 
sumed.; injection: Magnetized matter is continuously injected into 
the computational domain near the outer boundary 


the flow. Note that SPOl and PB03 performed 2D sim¬ 
ulations, while others did 3D simulations. Despite some 
differences in simulations, both of SPOl, INA03, and 
PB03 obtained a relatively steep density profile, p (x r~^. 
No obvious density peak is found in their simulations. 

HB02 obtained a similar density profile but with a re¬ 
markable peak at ^ 5 rg, which may be a remnant of 
the initial torus. Caution should be taken here, since 
this peak may be transient. Similarly, the density peak 
found at ~ 20 rg in KMS04, which mainly contributes 
to X-ray emissions, might be a direct consequence of the 
initial torus located at 40 Cg at the beginning of the sim¬ 
ulations. If this density peak is removed by long term 
numerical calculations, we expect that the large X-ray 
excess may disappear and the resultant spectra can be 
consistent with the observations, even if we do not arti¬ 
ficially restrict the emission region. 

To see if this is the case, we made the following test. 
We performed 2D MHD simulations, starting from the 
similar initial condition to that of KMS04 but with the 
initial torus being put a much larger radius, around 
> 100 Tg. After a long time interval, we have confirmed 
that the broad density peak, a remnant of the initial 
torus, disappears. The resultant density structure in this 
2D model (at t = 8200) is plotted by the solid curve in 
the upper panel in Figure 3 with that of KMS04. The 
density profile of our 2D model is roughly proportional to 
at r > 7 Tg and roughly agrees with the simulations 
by SPOl, INA03 and PB03. 

We, next, calculate the spectrum based on our 2D 
model and plot it in Figure 3 (lower panel). The adopted 
density parameter is po = 3.0 x 10“^®gcm“^, and the size 
of the emitting region is i? = 30 rg. The calculated elec¬ 
tron temperature is shown in Table ^ Clearly, our 2D 
MHD model can nicely fit the observed data points dur¬ 
ing the flaring state. Thus, it is very likely that if we 
continue MHD simulations from the data by KMS04 fur¬ 
ther, we may be able to fit the observational data. We 
need caution, however, because behavior of the magnetic 
fields in 2D and 3D simulations may be distinct. We 
cannot trust the 2D data so much. We need to perform 
longer 3D simulations to confirm our tentative conclu¬ 
sions. 

It is expected that the MHD models by INA03, SPOl, 
and PB03 can also reproduce the observed SED dur¬ 
ing the flaring state, since these models have similar 
density profiles to those of our 2D model (see the up¬ 
per panel of Figure 3). We demonstrate it by calculat¬ 
ing the spectrum from the regions at r < 30 rg of the 
spherical flow model, which is constructed based on the 
3D MHD simulations of INA03 (see the dotted curve in 
the lower panel of Figure 3). In this model the den¬ 
sity and proton temperature profiles are set to be p = 
7.0x 10“^^(r/rg)“^gcm“^ and Tp = 8.6x 10^^(r/rg)“^K, 
respectively. The electron temperature and the plasma- 
(3 are assumed to be constant in the calculating box, 
Te = 5 X 10®K and /3 = 3. The INA03 model shows 
a similar SED with that of our 2D model, and is con¬ 
sistent with the observed data in the flaring state. The 
other MHD models by HB02, MMMOl, and MM03, on 
the other hand, have flatter density profiles, p cx r““ 
with a = 0 — 1/2 (see Table |21l. Hence, they will produce 
an X-ray excess as was in the case of our simulation. 
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Fig. 3.— Upper panel: The normalized density profile given by 
our 2D MHD simulations (thick solid) and KMS04 (dashed). The 
density peak in KMS04 is around 20 rs, in contrast, it is near the 
black hole (r ~ 7 rs) in the case of our 2D simulations. The thin 
solid curves indicate the schematic density profile given by IN AOS 
(3D simulations), SPOl, and PB03 (2D simulations). They showed 
a similar density profile, p oc r“^, with our 2D MHD simulations 
at the outer part (r > 10 rg). Lower panel: The solid curve is 
the spectra of our 2D MHD model. We consider the emission from 
the extended region of r < 30 rg. The normalization of density 
is set to be po = 3.0 X 10“^®gcm“^, and the calculated electron 
temperature is shown in Table 0 The dotted curve represents 
the spectra of the INA03 model. Here, we use the spherical flow 
model, which is constructed based on the 3D MHD simulations of 
INA03. In this model, the density and proton temperature profiles 
are p = 7.0 X 10“'^® (r/rs)“^g cm“® and Tp = 8.6 X 10^^ (r/rs)“^K, 
respectively. We set Te = 5 X 10®K and plasma-/? = 3.0. Our 
2D MHD model as well as the INA03 model can reproduce the 
observed data during the flaring state. 

4.1.2. Moderate radiative cooling 

A distinct way of possibly producing a steeper density 
profile is to introduce moderate radiative cooling (Mi- 
neshige et al. 2002). By the terminology of RIAF it 
means that radiative loss is inefficient. However, it may 
happen that radiative cooling is only partly and tran¬ 
siently important. If we set the cooling timescale to be 
comparable to the accretion timescale, we expect that 
only dense parts can efficiently be cooled. Since the den¬ 
sity is generally higher in the innermost region, cooling 
results in further mass concentration towards the center, 
a preferable condition to fit the observations. However, 
this is only a possibility and we need 3D MHD simula¬ 
tions to confirm this idea. 

Another possibility is that the outer portions of the 
disk may become radiation efficient flow. Even if the 
density profile does not change, the emission from the 
outer parts can be reduced, if the electron temperature 
there were much lower than that of the inner part. In 
KMS04 and most of other MHD simulations, the one- 
temperature plasma is assumed and the cooling of pro¬ 


tons is neglected. However, the proton thermal energy 
may be converted into the electron energy due to the ef¬ 
fective Coulomb collisions and be radiated away in the 
outer zones. Further study requires 3D MHD simulations 
coupled with the electron energy equation and radiative 
cooling. 

4.2. Quiescent state 

The case of the quiescent spectrum will be much more 
difficult to reproduce by the MHD model. We have seen 
that the calculated spectral slope was much steeper than 
that of the X-ray observations in the quiescent state of 
Sgr A*. We can reduce the X-ray flux but cannot easily 
change the spectral slope in the X-ray range. We find 
two possibilities to make a flatter spectral slope. 

4.2.1. Nonthermal electrons 

This difficulty may be removed, if emission from the 
nonthermal electrons dominates over that from thermal 
ones. By considering emission from nonthermal elec¬ 
trons, Yuan, Quataert, & Narayan (2003) nicely repro¬ 
duced the observed spectrum in the quiescent state of 
Sgr A*. [However, they prescribe an arbitrary energy 
distribution of the nonthermal electrons with more free 
parameters in order to give a good fit to the observa¬ 
tions.] This is a very attractive possibility, but it is hard 
to prescribe an energy distribution of nonthermal elec¬ 
trons because of a lack of good theory of particle acceler¬ 
ation. Thus we need to introduce parameters describing 
the electron energy distribution and the results strongly 
depend on these parameters. 

4.2.2. Compton scattering in disk coronae 

The spectral slope at the X-ray band would be flatter if 
inverse Compton scattering is more enhanced, compared 
with other processes, and becomes dominant in the X- 
ray range. Note that the calculated steep spectral slopes 
shown in Figure 1 are due to bremsstrahlung emissions 
from the outer parts. The relative importance of the 
Comptonization inside the corona is proportional to the 
column density, ~ pH (with H being the half-thickness 
of the corona), while the bremsstrahlung emissivity is 
proportional to p^H. Thus, the presence of an extended, 
tenuous corona may resolve the problem. However, the 
formation mechanism of the corona is still an open issue 
and it is not easy to probe this possibility 

5. CONCLUSIONS 

In order to test the 3D MHD flow simulation as a 
model of optically-thin, high-temperature accretion flows 
through comparison with observations, we studied the 
spectral properties of magnetized accretion flows based 
on the 3D MHD simulation data. Surprisingly, such com¬ 
parative studies have not been well investigated so far in 
spite of a large number of MHD simulations having been 
performed recently. We summarize our results as follow¬ 
ings: 

1. We found that the MHD model cannot reproduce 
the observed SED in the flaring state of Sgr A* 
without substantial modification. If we use the 
data by KMS04, we need to require that the emis¬ 
sion region is restricted to be compact (r < 10 rs). 
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but this could be due to the particular initial con¬ 
dition. Some other MHD models with a steeper 
density profile {p oc r~^) and without a broad den¬ 
sity peak (r > lOrg) can fit the observations, if 
the emissions from regions of r < SOrg mainly con¬ 
tributes to the SEDs. 

2. The MHD flow can not generally reproduce the 
spectrum (spectral shape, in particular) in the qui¬ 
escent state. Significant contributions by nonther- 
mal electrons should resolve this issue. It is also 
possible that observational X-ray data may con¬ 
tain emissions from regions other than the vicinity 
of the black hole. 

3. The MHD flow predicts substantial and incoherent 
time variations in the emergent spectrum. The pre- 
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